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Abstract 

In a sample of 5.3 x 10 7 0-decays observed with the KLOE detector at the Frascati 
0-factory DA<J>NE we find 605 r]7t°'-f events with r\ — > 77 and 197 ry7r°7 events with 
r] — > 7r + 7r _ 7T°. The decay r/7r 7 is dominated by the process 0^ao7- From a fit to 
the 7?7T° mass spectrum we find BR(0^a o (98O)7)= (7.4 ± 0.7) x 10" 5 . 

Key words: e + e~ collisions, <p radiative decays, Scalar mesons 
PACS: 13.65.+i, 14.40.-n 



There is no clear understanding of the ao(980) and /o(980) mesons in the quark 
model. It has been argued that these mesons might not in fact be qq states but 
rather 4-quark states {qqqq) [1] or KK molecules [2]. The amplitude for the El 
transition 0— >ao7 to a state ao with J p (ao) = + is proportional to k x /(m 2 ), 
with k the momentum and m the mass of the ao- The differential decay width 
dr/dm for 0— >ao7 has the form k 2 x F(m 2 ), where F is a Lorentzian for the 
ao times a damping factor depending on the <fi and ao wave functions. The 
shape of the a signal in decay is therefore grossly distorted, acquiring a 
large tail at low m. 

According to different interpretations the — > a 7 branching ratio can range 
from 10~ 5 for qq and KK to 10~ 4 for qqqq [3]. The ratio BR(0^/ O 7) /BR(0^a O 7) 
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also depends on the structure of the scalars [4]. Chiral perturbation theory also 
attempts to estimate the amplitude and the distortion of the ao line shape ob- 
served in meson decays [5]. Production of the ao meson, followed by ao— >777r° 
dominates the final state rjTi ^ in 0-decays. A small contribution is due to 

We present in the following a study of the decay <f)^r]7r '-f performed with the 
KLOE detector [6] at the Frascati 0-factory [7] DA$NE. DA$NE delivered 
about 16 pb -1 during the year 2000, operating at a total energy of 1020 MeV. 
About 5.3 x 10 7 (p mesons were produced and collected by KLOE. We searched 
for two types of events. Type 1 events have five photons from 0— >?77r 7 with 
rj— >77_ Type 2 events have five photons plus two charged tracks corresponding 
to the same initial decay but with r]-^ii + ii^7i . Observation of type 1 events 
has been reported by SND [8] and CMD-2 [9] experiments at Novosibirsk. In 
this work we report the first observation of type 2 events. While the yield for 
type 2 events is lower, there is no physical background. 

The KLOE detector consists of a large cylindrical drift chamber surrounded by 
a lead-scintillating fiber electromagnetic calorimeter. A superconducting coil 
provides a 0.52 T field. The drift chamber [10], 4 m in diameter and 3.3 m long, 
has 12,582 all-stereo tungsten sense wires and 37,746 aluminum field wires. The 
chamber shell is made of carbon fiber-epoxy composite and the gas used is a 
90% helium, 10% isobutane mixture. These features maximize transparency 
to photons and reduce Kl — > Kg regeneration and multiple scattering. The 
position resolutions are ~ 150 /im in the coordinate transverse to the wire 
direction and « 2 mm in the longitudinal one. The momentum resolution 
is u(pi_)/p± ~ 0.4%. Vertices are reconstructed with a spatial resolution of 
~ 3 mm. The calorimeter [11] is divided into a barrel and two end-caps, for a 
total of 88 modules, and covers 98% of the solid angle. The modules are read 
out at both ends by photomultipliers; the readout granularity is ~4.4x4.4 cm 2 , 
for a total of 2,440 "cells." The arrival times of particles and the positions in 
three dimensions of the energy deposits are determined from the signals at the 
two ends. Cells close in time and space are grouped into a calorimeter cluster. 
The cluster energy E C l is the sum of the cell energies. The cluster time tch 
and position r*cL are energy weighed averages. Time and energy resolutions 
are a E /E = 5.7%/^E(GeV) and a t = 57 ps/y^GeV) © 50 ps, respectively; 
space resolution is ~ 1 cm in all the three coordinates. The detector trigger [12] 
uses calorimeter and chamber information. 

We first consider the case of <p — ► rj-K '-/ decays in which t] — > 77. These 
events are characterized by five prompt photons. A photon is detected as a 
calorimeter cluster with Eq^ > 3 MeV, with no track pointing to it, and 
satisfying \t CL — |r*cL|/c| < 5a t (E CL ). Photons with (cos^l > 0.93 are rejected 
to avoid machine background. Exactly five prompt photons are required. The 
main background processes are: 
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(1) — > 7r°7r°7 dominated by 0— >/o7; 

(2) e + e~ — > cutt with a; — * 7r°7 

(3) 0— >?77 with 77— >77 

(4) 0— >?77 with 77 — > 7r°7r 7r 

Process (3) can mimic five photon events due to energy cluster splitting or ac- 
cidental background while in process (4) photons can escape detection. Back- 
ground from (fr^KsKi with Kg— >ir ir and Kl interacting in the calorimeter 
is negligible after requiring the sum of the energy of the five prompt photons 
to be greater than 700 MeV. About 1.5xl0 4 events survive cuts. A first kine- 
matic fit, in which 4- momentum conservation and t—\f\/c = for each prompt 
photon are required, is performed. Background from processes (3) and (4) is 
reduced by requiring x 2 /ndf<3. Events with a 77 pair having an invariant 
mass close to the rj mass within 30 MeV are retained. Then events with one of 
the three remaining photons having energy greater than 340 MeV are rejected 
to reduce background (3); 2.5 x 10 3 events remain after cuts. 

The best photon pairing is found by matching the invariant masses of the 77 
pairs to the intermediate particles masses, either (i) 77 and ir° or (ii) 2 7r°. A 
second kinematic fit is then performed, with constraints on the masses of the 
intermediate particles, in both hypotheses. For hypothesis (i) we retain events 
with x 2 /ndf<3. The sample is still dominated by 7r°7r°7 events from processes 
(1) and (2). These two contributions can be discriminated from the signal by 
exploiting the result of the fit with hypothesis (ii). For each of the two 7r°7 
combinations, an invariant mass M nl is obtained; in Fig. 1 the absolute value 
of their difference AM^ is plotted versus the 7r°7r° invariant mass M n7T , both 
for data and Montecarlo (MC) events. The um° contribution is represented by 
the curved band, while — > 7r°7r°7 events are mainly located at high values 
of M nn . In order to select a clean — > r]7T 0/ ~f sample, the region below the 
solid curve of Fig. 1 and with M n7T < 760 MeV is retained. The final sample 
consists of 916 events. The spectrum of the rjir invariant mass M vn is shown 
in Fig. 2 together with the expected distribution for the background. In Fig. 
2 the distribution of cos# 7 of the unassociated photon is also shown, exhibit- 
ing the expected l+cos 2 # 7 behaviour. The efficiency for the identification of 
signal events is evaluated by applying the whole analysis chain to a sample of 
simulated — > rjir '-/ events generated with a uniform M V7T distribution. The 
selection efficiency as a function of M vn is shown in Fig. 3. The average over 
the whole mass spectrum is 32.4%. The background rejection factors and the 
expected number of events, as obtained by MC simulation, are given in Tab. 
1. The total number is 309±20 background events. 

We next consider the case of r)7T ~f in which rj — > 7r + 7r~7r°. This decay chain 
gives a final state with 2 charged pions and 5 prompt photons. This signature 
is unique among the possible final states so that the main background sources 
come from final states with similar topologies and much larger branching ra- 
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Fig. 1. Invariant mass difference between the two 7r°7 combinations versus the 
7r°7r° mass in hypothesis (ii) (see text); (a) data, (b) MC (j> — > r/7r°7, (c) MC 
e + e~ — > unr° — > 7r°7r°7, (d) MC <f> — > 7r 7r°7. The solid lines delimit the selected 
region. 




Fig. 2. i] — > 77 sample. Left: ??7r invariant mass spectrum (points), residual back- 
ground contribution (dashed histogram); Right: cos# 7 distribution of the unassoci- 
ated photon, comparison between data (points) and MC signal (dashed histogram) 
and MC background (dotted histogram). 



tios. The most significant backgrounds are: 



- — > 777 with 7] — > 7t + 7t~7t° (2 tracks and 3 photons); 

- e + e" — > wit with uj — > 7r + 7r~7r° (2 tracks and 4 photons); 

- — > KsKl with a prompt i\~z, decay and a combination of — > 7r + 7r~ and 
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(MeV) 

Fig. 3. Selection efficiency as a function of M„ n for type 1 and type 2 final states. 

Kl — > 7T 7r 7r° or Ks — > n°7r° and Kl — > ixlv or Kl — > 7r + 7r~7r° resulting in 
2 tracks and 4 or 6 photons. 

The signal events are selected by requiring a vertex close to the interaction re- 
gion with two tracks of opposite charge, and five prompt photons with E 1 >10 
MeV and |cos#| < 0.93. The surviving events (7.1 x 10 3 ) enter a two-step kine- 
matic fit procedure: a first fit (fitl) with energy /momentum conservation at 
the vertex, and a second fit (fit2) with tt° and rj mass constraints, where the 
combination resulting in the higher x 2 probability is selected. The invariant 
mass distributions for all 77 and 7r + 7r~77 combinations are shown in Fig. 4(a) 
and (b) for the events surviving fitl. A clear 77 signal already emerges at this 
stage over a large combinatorial background due to residual un° and K$K L 
events (see the peak close to the uj mass in Fig. 4(b)). After fit2 the residual 



Process 


Rejection factor 


Expected events 


e + e~ — » lott — > 7r°7r°7 


140 


54± 6 


<j> -» vr 7r 7 


40 


152±16 


(f> — > 777 — > 777 


6xl0 4 


5±2 


4> — > 7/7 — > 7r 7r°7r°7 


2.5xl0 3 


98±10 



Table 1 

Rejection factors for the background processes contributing to the five photon fi- 
nal state. The errors on the number of expected events include MC statistics and 
systematics from cross section uncertainties. 
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Fig. 4. Distributions of the invariant mass of all 77 and 7r + 7r~77 combinations for 
the events surviving fitl ((a) and (b)) and fit2 ((c) and (d)). The distributions are 
compared with the MC ones (histograms) for the signal only. 

background is reduced to 4±4 events from lv7t° and less than 8 events from the 
other possible background modes. Fig. 4(c) and (d) shows the 77 and 7i + 7i~ r yy 
invariant mass distributions at this stage compared with the MC distributions 
for the signal normalized to the same number of events. The comparison shows 
good agreement. The distribution of M vn for the 197 events found is shown in 
Fig. 5(a). The angular distribution of the radiated photon, Fig. 5(b), agrees 
with (l+cos 2 6» 7 ). 

The efficiency (see Fig. 3) as a function of M w is found by using MC sim- 
ulation. Both photon and track efficiencies are corrected using data control 
samples. (K$ — > ir + ir~ for tracking and e + e~7 for 7's). The overall efficiency, 
on average 19%, is dominated by the probability for at least one of the seven 
particles to go out of acceptance. The resolution in M V7T is 4 MeV over the en- 
tire mass range. A non Gaussian tail of about 10% is present, due to incorrect 
photon pairings. 

The — > 77"7r°7 branching ratio is obtained independently for each sample, 
normalizing the number of events after background subtraction, N — B, to the 
4> cross section a^, to the selection efficiency e, and to the integrated luminosity 
L: 

BRU -> mr°7) = ^LzA — x — \ — . (1) 

VV 1 U ex BR(r) z) (T^xL K J 



where % = 77, tt + ti tt° and BR(?7— n) is from ref. [13]. The luminosity is mea- 
sured using large angle Bhabha scattering events. The 0^ measurement is 
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Fig. 5. 77 — > 7r + ?T _ 7T sample, (a) Distribution of the 777r invariant mass for the final 
sample of events; (b) distribution of |cos(# 7 )| for the final sample of events. The 
dashed histogram is the MC, and the curve superimposed corresponds to a (1 + 
cos 2 (6> 7 )) dependence. 

obtained using the — > 777 — > 777 decay [14]. 



We obtain for the sample in which 77 — > 77 

BR(0 -> 777^7) = (8.51 ± 0.5W ± 0.57 syst ) x lO" 5 
and for the sample in which 77 — > 7r + 7r~7r° 



BR(0 -> r77r°7) = (7.96 ± 0.60 stat ± 0.40 



syst,, 



x 10" 



The two values are in agreement. Our results also agree with those from SND 
(8.8±1.4±0.9) x 10~ 5 [8] and CMD-2 (9.0±2.4±1.0) x 10~ 5 [9]. Contributions 
to the uncertainties are listed in Tab. 2. 

We estimate the contribution of 0— *-ao7 from a simultaneous fit of the two M m 
spectra. Two contributions are considered in the fit: (a) 0— >p°7r°, p°— ^777 and 
(b) 0— >a 7, ao^rjir . The contribution from the continuum process e + e~ — »■ 
ujtt , uj — > 777 is negligible due to the low cross section. The M vn spectrum for 
(a) is taken from VMD calculations [15,16]. For (b) we use the formulation of 
ref. [17] based on coupling to a charged kaon loop: 



d r(0 -> a 7 -> 777r°7) 2M% r r 0ao7 r ao7?7r o 



dM 



7J7T 



7T 



ID I 2 



(2) 
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Type 1 


Type 2 


statistics 


0.43 


0.58 


background subtraction 


0.28 


0.15 


efficiency 


0.51 


0.30 


BR(t] -> x) 


0.05 


0.14 


luminosity 


0.17 


0.16 


4> cross section 


0.17 


0.16 



Table 2 

Contributions to the uncertainties in BR(<p — ► rjtr "/) measurement (10~ 5 units). 

where r^ ao7 is related to the coupling g^ oK+K -/4:7r, T aor]7T o to the coupling 
9a 7]7r° I '^ 7r ' anc ^ D ao is the inverse a propagator including finite width correc- 
tions. The model assumes point-like couplings. 

The efficiency and resolution functions, including distortions from incorrect 
photon pairings are folded into the theoretical distribution. The experimental 
spectra after background subtraction (Fig. 6) are simultaneously fitted, setting 
M ao =984.8 MeV, from ref [13]. The free parameters of the fit are the branch- 
ing ratio for contribution (a) and the two coupling constants above. We find 
BR(0^p°7r°, p°-^7)=(0.5 ± 0.5) x 10~ 5 , g^ K + K -/^ = (0.40 ± 0.04) GeV 2 , 
and gaow/9a K+K- — 1-35 ± 0.09. The contribution from 0— >a 7 is dominant 
and that from p°7r° is consistent with zero, in agreement with VMD calcula- 
tions [16]. By integration we find: 

BR((p -> a 7, a -> r]n ) = (7.4 ± 0.7) x 10" 5 (3) 

The fit is good, with x 2 /ndf=27.2/25. Fig. 6 shows the fit and the ao contri- 
bution. Interference between (a) and (b) has been considered [15] giving no 
significant effect in the fit results. 

The result for g aoV n/ g ao K+K- can be compared with the value 1.05±0.06 given 
by experiment BNL852 [18] and with the Crystal Barrel results (1.07 or 0.93 
depending on the analysis [19]). It can also be compared with expectations 
based on the standard qq model (1.51 according to [17]) and on the qqqq model 
(0.93 according to [1]). 

Finally, the results are combined with those obtained in the analysis of <ft —>■ 
7r°7r°7 [20] done on the same data sample: 

- BR((p -> / 7 -> 7r°7r°7) = (1.49 ± 0.07) x 10" 4 ; 

- gj oK+K _/^ = 2.79 ±0.12 GeV 2 . 
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Fig. 6. Result of the combined fit: (a) comparison of data (exp. points) vs. fit 
(histogram) for (a) type 1 (b) type 2 samples and (c) plot of the theoretical curve 
for the do contribution with the parameters extracted from the fit. 



Multiplying the branching ratio above by a factor of 3 to account for fo — > 
7r+7r~ decays BR(0— >f j) is then (4.47 ± 0.21) x 10~ 4 , and the ratio of the 
two branching ratios is 



_BR(0_^/ot) 
^-BR(0^a o7 )- 6 - 1±0 - 6 

The ratio of the two couplings to the KK system is 



(4) 



R g2 = 4^ = 7.0 ± 0.7 

9a KK 



(5) 



These ratios are useful for shedding light on the possible sizes and mixing or 
other structure questions regarding these mesons [4,21]. 
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